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Pollen per cone 


Sonora 
Dallas 
San Marcos 
Junction 


472.000 pollen grains/cone 

402.000 pollen grains/cone 

374.000 pollen grains/cone 

363.000 pollen grains/cone 



Cone and pollen production for 
representative trees 



Cones/tree 

Total pollen 
potential 

Santa Fe - LCP 

52,808 

1.53 x 10 10 

Santa Fe - HCP 

646,496 

1.87 x 10 11 

Jemez Springs - HCP 

269,946 

7.83 x 10 10 


Summary Statistics for 2010 and 2011 Juniperus 
pinchotii pollen seasons at source 


Location 

Average daily 
concentration 
Pollen 
grains/m3 

Peak daily 
concentration 
Pollen 
grains/m3 

Date of 
peak 

Peak hourly 
concentration 
Pollen grains/m3 

Time of 
peak hour 

Date of 
peak hour 

2010 

Erick, OK 

337 

5,563 

25-Oct 

15,898 

10:00 AM 

25-Oct 

Sonora, TX 

286 

3,019 

25-Oct 

12,152 

10:00 AM 

25-Oct 

San Angelo, TX 

653 

5,542 

1-Nov 

10,195 

Noon 

1-Nov 

2011 

Erick, OK 

12 

214 

16-Oct 

800 

10:00 AM 

16-Oct 

Sonora, TX 

64 

428 

26-Oct 

2,422 

10:00 AM 

31 -Oct 

San Angelo, TX 

60 

493 

26-Oct 

2,747 

4:00 AM 

26-Oct 

Quanah, TX 

22 

190 

19-Oct 

703 

8:00 PM 

19-Oct 


Influence of preseason meteorological 
variables for Juniperus ashei 

► Start date 

Significantly correlated with mean monthly temperature 
in December (r = 0.467, p = 0.038) and November 
rainfall (r = 0.468, p = 0.038) 

► Cumulative Season Total (CST) 

► Significantly correlated with mean maximum temperature 
in December (r = 0.4740, p = 0.035) 


Correlation of average daily Juniperus pollen 
concentration with meteorological variables from 1987 

to 2006 


Meteorological Variable 

Main season r 

Max daily temperature 

0.607*** 

Min daily temperature 

0.371*** 

Mean daily temperature 

0.546*** 

Rainfall 

-0.143*** 

Rainfall (1 day lag) 

-0. 164*** 

Rainfall (2 day lag) 

-0.069 

RH 

-0.282*** 

Mean wind speed 

0.117** 

Sunshine 

0.257*** 


* p <0.05 ** p<0.01 *** p<0.001 


Multiple regression model for main season pollen 
concentration.and meteorological variables 


R 2 

Meteorological variable 

Beta value 

P 

0.379 



<0.001 


Max daily temperature 

0.583 

<0.001 


Rainfall 

-0.071 

0.014 


Rainfall (1 day lag) 

-0.064 

0.027 


Mean wind speed 

0.033 

0.253 


Multiple Regression Model for Start Date and 
Preseason Meteorological Variables 


R 2 

Meteorological variable 

Beta value 

P 

0.736 



<0.003 


November rain 

0.795 

0.001 


December mean T 

0.566 

0.004 


J anuary rain 

0.559 

0.009 


October rain 

-0.500 

0.016 


November mean T 

0.502 

0.032 


August rain 

0.316 

0.104 


Multiple Regression Model for CST and 
Preseason Meteorological Variables 


R 2 

Meteorological variable 

Beta value 

P 

0.639 



<0.019 


December maximum T 

0.754 

0.002 


J une rain 

0.469 

0.021 


December rain 

-0.384 

0.045 


J anuary mean T 

-0.237 

0.211 


September rain 

-0.285 

0.148 


J anuary rain 

-0.208 

0.262 


Red Cedar Encroachment 


> Oklahoma has 17 million acres of prairie, shrub land, 
cross timbers forests and other forests 

> 1950: 1.5 million acres with cedar problems 

> 1985: 3.5 million acres with cedar problems 

> 1994: 6.3 million acres with 50 trees /acre and 2.5 milli on 
acres with 250 trees/acre - 37% loss of native ecosystems 

> 2001: 8.0 million acres with 50 trees /acre and 5.0 milli on 
acres with 250 trees/ acre - this represented a 47% loss of 
native ecosystems 

> 2013 projection: 12.6 million acres with 50 trees/acre and 
8.00 million with 250 trees/ acre 


Meters AG L Source* at multiple locations 


NCJACTWSFWnODEL 

Forward trajectories starting at 1800 UTC 08 Jan 1 1 

12 UTC 08 Jan NAM Forecast Initialization 
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Job ID ll£76 Job Start Sa: Jan 0 14 48 23 UTC 201 1 

Source 1 Jal PBS Ion -101 50 heigh I IWrtlAGl 

Trajectory Direction Forward Duration 43 hrs 
Vertical Moll on Calculation Method Model Vertical Velocity 

Meteorology 12GGZ 8 Jan 2011 NAM 12 km 


Northerly Flow 


NOAA HYSPLIT MODEL 

Forward trajectories starting at 1800 UTC 29 Dec 10 

00 UTC 29 Dec NAM Forecast Initialization 
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Job ID 13972 Job Start Wed Dec 29 06 3 7 24 UTC 2010 

Source 1 let 29 5 Ion -10150 height lOOmAGl 

Trajoctory Direction Forward Duration 48 hrs 

Vertical Motion Calculation Method Model Vertical Velocity 

Meteorology 0000Z 29 Dec 2010 NAM 12 km 


NOAA HYSPLIT MODEL 
Forward trajectories starting at 1800 UTC 02 Jan 1 1 
GDAS Meteorological Data 
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Job Start: Frl Feb 18 05:37:23 UTC 201 1 
- 1 01 .50 hetghl : 100 m AGL 

Trajectory Direction: Forward Duration: 48 hrs 

Vortical Motion Calculation Method: Model Vertical Vetodty 

Meteorology: OOQOZOI Jan 2011 - GDASl 


Meters AGL Sourced at multiple locations 


NOAA HYSPLIT MODEL 

Forward trajectories starting at 1 800 UTC 02 Feb 1 1 
00UTC02Feb NAM Forecast Initialization 
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Job ID: 15395 Job Start: Wed Feb 2 03:25:24 UTC 2011 

Source 1 lat.:29.5 Ion.: -101.5 height: 100 m AGL 

Trajectory Direction : Forward Duration : 48 hrs 

Vertical Motion Calculation Method: Model Vertical Velocity 

Meteorology: 0000Z 2 Feb 2011 - NAM 12 km 


Southerly Flow 


NOAA HYSPLIT MODEL 

Forward trajectories starting at 1 800 UTC 01 Feb 1 1 
00 UTC 01 Feb NAM Forecast Initialization 
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Job ID: 13218 Job Start: Tue Feb 1 04:49:59 UTC 201 1 

Source 1 lat. : 29.5 Ion.: -101.50 height: 100 m AGL 


Duration : 48 hrs 
hod: Model ' ‘ 

Meteorology: 0000Z 1 Feb 2011 - NAM 12 km 


Trajectory Direction: Forward 

Vertical Motion Calculation Method: Model Vertical Velocity 


NOAA HYSPLIT MODEL 

Forward trajectories starting at 1800 UTC 10 Jan 1 1 

00 UTC 10 Jan NAM Forecast Initialization 
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Job 10: 13446 Job Siart Mon Jan 1 0 04 i6:l£ UTC £011 

Source 1 lal 20 5 Ion -tt1» heigh t 1 n m AGL 

Trajectory DinMiOrt Forward Duration 48 hrs 

Vertical Motion Calculation Method Model Vertical Velocity 

Mcieorology 03Q0Z 10 Jan £&11 NAM 12 km 



Meters AGL Source* at multiple locations 




Incursion of J. ashei Pollen into Tulsa 

► Each winter since 1980 J. ashei pollen has been registered 
by our Tulsa air samplers 

► Levetin and Buck, Annals of Allergy, 1986. 

► Levetin, Aerobiologia, 1998 

► Rogers and Levetin, Int J Biometeorol, 1998 

► Van deWater and Levetin, Grana, 200 1 

► Van deWater et al, Int J Biometeorol, 2003 

► Pollen recorded on 48% of the days in Dec and Jan (range 
I I % to 79%) 

► Concentrations typically low; however, 'Very High ” 
concentrations have been registered on several occasions 
(based on National Allergy Bureau level of very high 

► 1 500 pollen grains/m 3 ) 


25 Oct 2010 incursion of Jurdperus pollen into the 
Tulsa atmosphere with a average daily concentration 
of 106 pollen/ m 3 



Trajectory analysis indicates the pollen originated in 
southwest Texas approximately 8 hours earlier 




NCAA HYSPLI I MODLL 

Forward trajectories starting at 1600 UTC 25 Oct 1 0 

GDAS Meteorological Data 


This is not a NOAA product. It was produced by a web user. 

Job ID: 349594 Job Start: Thu Feb 1 6 20:06:47 UTC 201 2 

Source 1 lat.: 30.26 Ion.: -1 01 .53 height: 500 m AG L 
Trajectory Direction: Forward Duration: 24 hrs 
Vertical Motion Calculation Method: Model Vertical Velocity 

Meteorology: QQQQZ 22 Oct 2010 - GDAS1 
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Juniperus pinchotii concentrations at Sonora 
and San Angelo confirm the trajectory model 

Hourly concentrations of Juniperus pinchotii pollen 
25 Oct 2010 - Sonora, TX 
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Hourly concentrations of Juniperus pinchotii pollen 
25 Oct 2010 at San Angelo, TX 
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Estimated weight per pollen grain (J. monosperma) after 2 hrs across 
the range of relative humidity levels at 20° C 


While J. monosperma and J. pinchotii were 
larger in size, their percent weight gain was 
similar to that of J. ashei. 



Relative humidity 


Estimated weight per pollen grain (J. pinchotii) after 2 hrs across the 
range of relative humidity levels at 20° C 


Hygroscopic Weight Gain of Juniperus 

Pollen 


• Juniperus oshei, J. monosperma, 
and J. pinchotii were dusted on 
greased slides and exposed to a 
range of temperatures and 
humidity levels and weighed at 2 
hrs and 6 hrs. 

•7. ashei was weighed at 
temperatures 4° C, 15° C, and 
20° C at relative humidity levels; 
20% 40% 50% 75% 86% 97%. 


•Weight was not significantly 
affected by temperature or time. 


2.90E-06 


♦ 



Estimated weight per pollen grain (J. oshei) after 2 hrs across the 
range of relative humidity levels at 20° C 
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Pollen release potential Source Map/Mask 


(PRPSM_of_Jj) of a Juniper species “f is calculated as: 

PRPSM_of_Ji = Tj x Mj x Hj x Cj x Pj 


Where 

Tj = Number of Jj trees 
M, = Male/Female ratio of Jj 
Hj = HCP_LCP/AII ratio for Jj 
Cj = Cones per Jj tree 
Pj = Pollens per cone for Jj 

The number of trees of a Juniper species “i” per grid cell is calculated as 

Ti = GAP; x MODIS x TC 

Where 

GAPj = Fraction of Jj at 1 km grid (range 0-1) 

MODIS = MODIS derived percent tree cover per 1 km 2 grid cell (in fraction, range 0-1) 
TC = Tree count or number of trees. 


Ground truth (transect data) 


(a) Male to Female ratio 

(b) HCP.LCP to All ratio 

0 - Only enough cones to determine gender 

1 - Low Cone Production (LCP) tree 

2 - High Cone Production (HCP) tree 

(c) Cones per tree 

(d) Pollens per cone 
(d) Age 

(Height & edge effect) 


Juniper Species and Pollination Season 


Juniper Ashei (JJ is mostly found to be distributed over Texas and 
Oklahoma and pollinates during December to January. Thus, the 
dispersion of juniper pollens during December-January is mostly 
restricted to J a type.. 


Juiper Pinchotti (J p ) is mostly distributed over Texas and pollinated 
during late September-November. Thus, the dispersion of juniper 
pollens during this period is mostly restricted to J p type. 


Juniper monosperma (J m ) and Juniper scopulurum (J s ) are 

prevalent in New Mexico and pollinates during March-May period. 
Thus, the dispersion of juniper pollens during this period is mostly 
restricted to J m and J s type. 


Field data 

• Juniper Ashei 

• Juniper Pinchotii 

• Juniper monosperma and scopulorum 

Information Needed (to update mask): 

For all sampling sites: 

• HCP/LCP/0 

• Male/Female ratio 

• Number of trees (tree density) 

• Number of cones 

• Pollens/cone 


Pollen Source Mask (PREAM) 








Juniper Ashei 

December to January 


Juniper Pinchotii 

late September-November 


Juniper 

Monosperma & Scopulorum 

March -May 


Colorado Kansas 


Missouri 


Oklahoma 


New Mexico 


Arkansas 


Colorado 


Missouri 


Arkansas 


Louisiana 


Louisiana 
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Missouri 


New 


Arkansas 


Pollen Count at Source 


GAP derived Juniper distribution 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


GAP derived Juniper distribution 


Pollen Count at Source 


Pollen Count at Source 


GAP derived Juniper distribution 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


Colorado 


Kansas 


Missouri 


Oklahoma 


New Mexico 


Arkansas 


Louisiana 


Nebraska 


Nebraska 


Kansas 


California 


Texas 


California 


Kansas 


Texas 


Spatial resolution: ~1 km (990 m) 


PREAM 


Pollen Plume Simulation for Juniper Emissions 
For the period 15 December 2009 - 1 January 2010 

Run by Slobodan Nickovic, September 2012 


Atiriosphere Model Setup 

Model horizontal domain : Southwest US 

Model resolution : ~40 km 

Simulation period : 

15 December 2009 - 1 January 2010 

Boundary conditions : 1 degree global forecasts used to refresh 

• initial conditions every 24 hours 

• boundary conditions every 6 hours 





PREAM 

Pollen Plume Simulation for Juniper Emissions 
For the period 15 December 2009 - 1 January 2010 


PREAM_ ~ Pollen Re£ional_A tmosgheric Model 
Derived from DREAM (dust), modified to simulate pollen 

- 4 particles bins 

- PREAM is online driven by the NCEP/ETA 
Emission: 

- Viscous-sublayer parameterization 

- Emission dependent on friction velocity 
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PREAM: Surface pollen concentration (#grains/m3) 

Forecast base time: 00Z15DEC2009 valid time: 06Z15DEC2009 (+06) 
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Model Validation ( Juniper pollen count/ m 2 ) 
Observed Forecast 


San_Marcos,TX Dec-Jan 2009/2010 
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Status of Transitioning 
Pollen Data Into NM EPHT 

Progress in year 1 
(green oval): 

• Prepare interface 
for health client 
server 

• Prepare server for 
pollen data output 

Test server functions: 
(yellow oval): 

• Dependent upon 
receiving sample 
data from modeling 
team 

Activities for out 
years: (blue boxes) 



A 


Ingest and 
Process Test 
Data 


Prepare, Modify, 
and Test Client 
Server and Web- 
based Services 



Model Output 
Stored on AZ 


j 


t 

Develop 


Server 


DREAM NMM / 

PREAM Model 

J 


Interface for 
Health Client 
Server 



Remote Sensing Data 




Meteorological 
Forecast Data 







Post-processing Workflow 


Acquire Pollen Forecast Output Files (.txt) 


I Dream-pollen fctO6.03 .2006. 

|conc*l.e+18 time 6 12 18 24 

45.00 -127.00 -9999.000 -9999.000 -9999.000 -9999.000 

45.00 -126.75 -9999.000 -9999.000 -9999.000 -9999.000 

45.00 - 126.50 -9999.000 -9999.000 -9999.000 -9999.000 

45.00 -126.25 -9999.000 -9999.000 -9999.000 -9999.000 

45.00 - 126.00 -9999.000 -9999.000 -9999.000 -9999.000 

45.00 - 125.75 -9999.000 -9999.000 -9999.000 -9999.000 


44.25 -121.75 
44.25 -121.50 
44.25 -121.25 
44.25 -121.00 
44.25 -120.75 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


PREAM Post-Processing Script 


1 . Read text file contents into Python lists while 
parsing hourly columns into separate data 
structures. 

2. Convert Python lists into Python numPy 
Arrays. 

3. Reconfigure numPy array into correct array 
dimensions. 

4. Apply multiplication factor to each grid value. 
Export to GeoTiff using GDAL libraries. 

6. Write records to GSToRE geospatial database 
and enable WMS,WCS services. 


Store GeoTiffs on File System 


WMS/WCS Services 
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GSToRE Spatial Database 
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Newport, Rhode Island Sep 18-20, 2012 
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Pollen 4 Day Forecast - 201 0-1 2-29 00:00:00ZUTC 


Sample Animation for Juniperus ashei 
(Dec 29, 2010 to Jan 1, 2011) 


Newport, Rhode Island Sep 18-20, 2012 
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Report a Syndrome 


Demographics 


Syndrome information 


Clinical Findings: Chronic Lung Disease Exacerbation 

-Symptoms (Reported by Patient} 




Productive Cough? Yes No Nasal Discharge? Yes Nc 
Sore Throat? Yes ■ No Wheezing? i ~ j Yes Nt 

Underlying Lung Disease (Asthma/ CQPD)? Yes ' No 


Clinical Signs (from Physical Examination) 


4 


Temp(C) < 37,0 37.0 - 37.9 38.0 - 38.9 39.0 - 39.9 

Predominant Lung Findings Rales Wheezing Bilateral 

Skin Rash? Yes No Oral Lesions? Yes No 
Lymphadenopathy? Yes No Diffuse Localized 

Splenomegaly? Yes No Hepatomegaly? Yes No 


Laboratory and X-Ray Data- 




WBC Count < 3,000 S.000 - 10,000 10,001 - 15,000 > 15,000 

Platelet Ct < 50,000 50,000 - 100,000 100,001 - 150,000 >15 

Chest X-Ray: Normal Abnormal 

|_ Infiltrate |_| Hyperinflation |_J Cardiomegaly |_| Effusion 

02 Sat, (Room Air) Normal Abnormal 


Help 


Cancel 


Submit Report 
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Vkw Ail Species 


Join Us! 


We are looking for volu nteers to 
help us monitor plant and animal 
species found across the United 
States. Click 'Observe' to Join us( 




USA National Phenology Network 

The USA National Phenology Network brings together citizen scientists, government agencies, non-profit 
groups, educators and students of all ages to monitor the impacts of climate change on plants and animals in 
the United States. The network harnesses the power of people and the Internet to collect and share 
information, providing researchers with far more data than they could collect alone. 


Learn more about us 


What is phenology? 


Phenology refers to recurring plant and animal life cycle stages, or phenophasest* ? such as leafing and 
flowering, maturation of agricultural plants, emergence of insects, and migration of birds. Many of these events 
are sensitive to climatic variation and change, and are simple to observe and record. As an USA-NPN observer, 
you can help scientists identify and understand environmental trends so we can better adapt to climate change. 

Why is phenology important? 


USA-NPN News 


Phenology Feed 


Join the Conversation 


Phenoclimatology Position at UA 
Introducing the USA-NPN Video / 

Nature's Notebook: “How to Observe" Handbook and Training Videos-^ 
Phenology Special Issue in the Philosophical Transactions of the Royal Society 

USA-NPN Reports (including Strategic Plan and 2009 Annual Report} 

Call for Papers: 4th Annual PROSE in Tucson, AZ. October 2010 


► Recent Media Reports 

► Newsletter Archive 


Are you... ? 

New to phenology? 
Ready to start observing? 
One of our partners? 
Interested In creating a partnership? 

An educator? 
Interested In finding data to use? 

A media outlet? 


About Ue | Contact Us | Partners, | Privacy | Accessibility | Glossary | links j Madia Center 
USA Mationel Phenology hletwork, Matkmel Coord ineting Office 
1C5S E. Sccth St.. Tucaon, AZ SS721 
Admin Login 





